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Abstract 

Using generating functions techniques we develop a relation between the Hausdorff 
and spectral dimension of trees with a unique infinite spine. Furthermore, it is shown 
that if the outgrowths along the spine are independent and identically distributed, 
then both the Hausdorff and spectral dimension can easily be determined from the 
probability generating function of the random variable describing the size of the out- 
growths at a given vertex, provided that the probability of the height of the outgrowths 
exceeding n falls off as the inverse of n. We apply this new method to both critical 
non-generic trees and the attachment and grafting model, which is a special case of the 
vertex splitting model, resulting in a simplified proof for the values of the Hausdorff 
and spectral dimension for the former and novel results for the latter. 
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1 Introduction 



Random trees and their properties have been studied in several branches of probability the- 
ory, mathematical physics and science in general. Their applications span from phylogenetic 
trees pfl [2] and random folding of RNA molecules [3] to models of quantum gravity jU EJ El [7] , 
amongst others. 

In general terms a random tree ensemble is a set of trees together with a probability 
measure associated to it. We focus on a class of infinite, rooted trees which have the property 
that there is a unique, non-backtracking path from the root to infinity. This path will be 
referred to as the spine. In this case the ensemble is characterized by the distribution of 
the finite outgrowths from the spine. Such ensembles can arise in very different contexts: 
One example considered in this article is the equilibrium statistical mechanics model of 
simply generated trees, first introduced in [S], which describes trees with a local action. 
These random tree ensembles are related to critical and sub-critial Galton- Watson processes 
[9]. The second example considered is the attachment and grafting model [10] which is a 
specific case of the vertex splitting model [11] . This random tree ensemble is very different in 
nature from the simply generated trees in the sense that it arises from a growth process with 
non-local action which results in a non-equilibrium statistical mechanics model. An effect 
similar to the emergence of a unique infinite spine is also observed in triangulation models 
in quantum gravity where a unique large "universe" emerges with finite baby-universes 
attached, see e.g. [T2] . 

A simple observable of a random tree ensemble is its Hausdorff or fractal dimension 
defined as dh provided that the number of edges within a distance R from the root, denoted 
by | -Br |, scales as 

\B R \ « R dh as i? ^ oo. (1) 

The meaning of "~ " is that both sides are asymptotically the same in a sense which will be 
made precise in Section [2J Another notion of dimensionality of the random tree ensemble is 
the so-called spectral dimension. Let p(t) be the probability that a simple random walk on 
a tree which starts at the root at time is back at the root at time t. If 

p(t) « t' ds/2 as t ->■ oo (2) 

we say that the spectral dimension of the tree is d s . The spectral dimension of various classes 
of graphs has been studied by probabilists [121 [EH H51 EES] and by physicists, especially in 
the quantum gravity literature, see e.g. [I7J [111 [Ell [20] . Coulhon [211 Theorem 7.7] derived 
that one has 

7~~~~r <d s <d h (3) 
1 + 4 

for fixed graphs under certain regularity assumptions. As we will see in the following, random 
tree ensembles with an infinite spine are tight with the left-hand-side of the inequality and 
we show below how one can obtain an improved inequality for those tree ensembles which is 
tight on both sides for many examples. 
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We develop methods which relate the diffusion of the random walk to properties of the 
volume distribution of the finite outgrowths. An important concept is the hull dimension 
defined as dh, provided that the total volume of the hull Br, the outgrowths from vertices 
on the spine within distance R from the root, scales as 

\8 R \ « R dh as R -»■ oo (4) 

see Figured) 

Using simple generating function techniques we begin by proving in Theorem [I] that the 
following bounds hold for trees with a unique infinite spine 

2 ' !h < d, < (5) 



1 + 4 l + d h 

which, in many examples of interest, tightens the general bound fl3]) for fixed graphs. In- 
tuitively, if the probability that the outgrowths reach height n decays fast enough with n, 
dh and d h should be close in which case (jSJ) gives a tight bound on d s . We make this a 
precise statement in Theorem [5] in the case when the outgrowths from different vertices on 
the spine are independent and identically distributed (i.i.d.). There, we prove that if the 
random variable \A l \, denoting the size of the outgrowths from a vertex on the spine, has 
a probability generating function (z^ A ^) = 1 — (1 — z) a l{l — z), where l(x) varies slowly at 
zero, then almost surely 

d h <-, d s < — . (6) 
a 1 + a 

Furthermore, we show that if the probability that the height of the outgrowths from a vertex 
on the spine exceeds n falls off as the inverse of n then one has equality in (J6j). 

Up to that point, our results are very general and we devote the rest of the article to 
applying the results to specific models. In Section [3] we calculate the spectral dimension 
of the so-called non-generic, critical phase of simply generated trees, giving an alternative 
proof to the recent results by Croydon and Kumagai [15]. We note here that our result 
is weaker in the sense that we define d s through the singular behaviour of the generating 
function of p(t), rather than from the asymptotic behaviour of p(t). In SectionHJwe calculate 
the Hausdorff and spectral dimension of the attachment and grafting model for a certain 
parameter range, extending results in [10J. We leave the generalisation to the full range of 
parameters as an open problem. 



4 



2 Hausdorff and spectral dimension of trees with infi- 
nite spine 



2.1 Trees with infinite spine and their ball and hull dimension 

Let us denote by T the set of all planar rooted trees, where the root has valency one. The 
planarity condition simply means that edges containing the same vertex are ordered around 
that vertex. One has T = T' U r°°, where T' is the set of finite such trees and r°° the set 
of infinite such trees. We often denote elements of T' by T, 7\ etc., while elements of r°° 
are called r, t% etc. For T G T', let |T| be the size of the tree, i.e. the number of edges in 
T. Let r denote the root of the tree. We consider in this article trees which have a unique, 
non-backtracking path from the root to infinity and call this path the spine. Denote the 
set of such trees by r^? C r°°. Furthermore, denote the vertices on the spine, ordered away 
from the root r, by Si, S2, S3, .... 

For a given r G T we denote by Br C t the ball of radius R around the root, this is the 
subgraph of r spanned by the vertices which have a graph distance from the root which is 
less than or equal to R. Another important concept in the article is the notion of the hull 
Br. For a r G V the hull Br denotes the subgraph which is the union of the spine from r up 
to vertex sr and all finite trees attached to it. The concepts of ball and hull are illustrated 
in Figure [U it is easy to see that Br C Br. 

A random tree ensemble (r, v) is a set of trees T equipped with a probability measure 
v. An important notion of dimensionality of a random tree ensemble, as discussed in the 
introduction, is the so-called Hausdorff or fractal dimension, which describes the growth of 
the size of a ball of radius R and can formally be defined as 

d h = hm — — . (7) 

R^oo hgR 

Throughout the paper we consider the following strong criterion for existence of dh- There 
exists an Rq > such that for R> i? 

R dh L(R) < \B R \ < R dh L(R), (8) 

where L(R),L(R) are functions which vary slowly at infinity. A real-valued, positive, mea- 
surable function L is said to vary slowly at infinity if for each A > 

lim ^ = 1. (9) 

We refer to [25] for some properties of slowly varying functions and in particular the fact 
that (ED implies (EJ). 

If for a random tree ensemble (r, v) one has that ((71) holds for v- almost all r G T then 
we call dh the quenched Hausdorff' dimension. Furthermore, one can also define the annealed 
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Figure 1: An illustration of the ball and the hull in a tree with a unique infinite spine 
started at the root r and labeled by vertices s±, S2, ... away from the root. The ball Br is the 
subgraph of all vertices within distance R from the root. The hull Br is composed of the 
spine up to distance R from the root and the finite outgrowths from its vertices sx, ...,sr. 
The finite outgrowths from vertex Sj are denoted by A 1 . 



Hausdorff dimension through the expected size of the ball of radius R 

\og(\B R \) v 



d 



H 



lim 

R-kx> log R 



(10) 



Furthermore, we say that the quenched Hausdorff dimension dh exists in the sense of (jHD if 
the latter is fulfilled for //-almost all r £ T. 

In analogy, we also define the quenched hull dimension, if for //-almost all r £ T 

log | Bjil 



dh = lim 



R->oo log R 



where existence is again provided through ((§]) analogous to the Hausdorff dimension. From 
the fact that Br C Br it is then also clear that one has dh < dh. 



2.2 Random walks on trees and spectral dimension 

In the following, we consider simple, discrete time random walks on trees i.e. walks which, in 
each discrete time step, jump to a nearest neighbour with uniform probability. For a given 
t £ r we denote by Q T the set of finite random walks on r. For a random walk u £ Q T we 
denote by \u\ the length of the walk and for t < \u\, u(t) denotes the vertex where the u is 
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located after t steps. We now formulate generating functions for the return probabilities of 
random walks on a given tree r. 

An important quantity is the generating function for the probability of first return to the 
root r given by 

P T (x)=Y^ P(Mt) =r\u>(0) =r,u(t') ^r,0<t' <t}){l-x)^ (12) 

t,u>£Q r 

and for all returns to a starting point which is given by 

Q T (x)= P (M*) =r\u(0) =r})(l-x)K (13) 

By decomposing a walk which returns to the root into a first return, second return etc., the 
two generating functions can be related through 

oo 1 

Q T {X) = Y,{PMT = — -J-. (14) 
n=0 T ^ ' 

The generating function Q T encodes all information of the asymptotic of the return 
probability and one can extract the spectral dimension from it. In particular, for a given 
t G T we define the spectral dimension through 

, s = 2 (l + liml2£^M), (15) 
\ x->a logx J 

provided that Q T (x) diverges as x — > 0. As for the Hausdorff dimension we also consider the 
strong criterion for existence, that there exists an xq G (0, 1) such that for x < Xq 

x- 1+ds/2 l_(x) < Q T (x) < x' 1+ds/2 l(x) (16) 

where l(x) and l(x) are functions which vary slowly at x — > + . Notice that if d s exists in 
the sense of (T161) . then the definition of d s as given in (Tl5il is equivalent to the more common 
definition 

i, = -2 lim bgK{"W = r| M (0) = r}) 

t^rOO log t 

under the assumption that P ({u)(t) = r \ u(0) = r}) is eventually monotonic in t. 

For an ensemble (r, v) we say it has quenched spectral dimension d s if (115j) is fulfilled for 
i/-almost all trees r G T. Furthermore, we define the annealed spectral dimension ds through 

d s = 2 [ 1 + lim W&W" 



log X 

provided that (Q T (x)) v diverges as x — > 0. 

To illustrate the concepts introduced in this subsection consider the following: 



7 



Example The easiest example is the (non-random) tree consisting of a single infinite spine 
which we denote by S. One has trivially that \Br\ = \Br\ = R and thus dh — dh — 1. To 
calculate the first return probability generating function Ps{x) note that the random walk 
leaves the root with probability 1, then at vertex si the walker returns to the root with 
probability 1/2 or leaves to the right with probability 1/2 and makes an excursion until its 
first return to Si. It then returns to the root with probability 1/2 or does another excursion 
with probability 1/2, and so on. Hence, one has for the generating function that 



From (1191) it follows that 
and thus from (fT^l) one has 



P S (x) = l-y/Z (20) 



Q s (x) = 4=- (21) 



Hence, using ( ITS]) one finds that d a = 1 as expected. 



2.3 Relating spectral and Hausdorff dimension 

The main result of this section is the following theorem which we prove in the two subsections 
below. 

Theorem 1 For any tree r G with a unique infinite spine for which dh, dh and d s exist, 
the following holds 

Idh , 2c?;, , , 

<d s < - T . (22) 



1 + d h ~ l + d h 

In the proof of this theorem, in the following two subsections, we use results developed 
in [22] in the context of generic trees and show that they can be employed in a more general 
context by introducing the hull dimension. When assuming existence of the dimensions we 
will in practice use the stronger criteria ([8]) and (fT6|) although the results also easily follow 
from assuming existence in the usual sence and (TTo^) . 

2.3.1 Upper bound on the spectral dimension 

For any tree r G T one can easily derive the following recursion relation [22] analogous to 
( |T9l) as presented in the example above 

PJx) = , (23) 
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where the n, Tu T denote the k T trees meeting the edge from the root to the unique vertex 
next to it. In [22] the following simple lemma is proven which follows straightforwardly from 
fl23l) by induction 

Lemma 2 For all finite trees T e T' one has 

P T {x) > 1 - \T\x. (24) 

To prove the upper bound on the spectral dimension we use the following lemma which is a 
reformulation of Lemma 7 in [22] 

Lemma 3 For all r G Vc? and R > 1 one has 

P T (x) > l-^-2xR-x\B R \. (25) 

Proof The proof follows directly the proof of Lemma 7 in [22] • Denote by P^(x) the 
generating function for the first return probability of walks which do not visit the vertex 
sr + \ on the spine. An induction proof using the recursion relation fT23|) then easily shows 
that 

1 - 

P T R (x) > l---xR-J2(l~ PA*)) . (26) 

i=l 

where A 1 denotes the union of the vertex Sj and the finite trees attached to it. Noting that 
P T {x) > P 7 R (x), using Lemma [2] and observing that \Br\ = ^2f =0 \A l \ + R completes the 
proof. 



To prove the upper bound in the theorem we note that from the existence of the hull 
dimension one has that \B T (R)\ < R dh L(R), where L(R) varies slowly at R — > oo. Thus one 
has from the proceeding lemma that 

QJx) > = (27) 

V 1 ; " 1 + 2xR 2 + xR^ +1 L(R) K ' 

Choosing R = \x^ l ^ 1+dhS> \ one has 

2.-1/(1+4) H x \ 

Qt{%) ~ l + /(a;) + 2x(4-i")/(i+4)/( 2; )' (28) 

where l(x) = 1/L( [ x - l ^ 1+ ^\ ) varies slowly at x — > + . Provided that d s exists and from 
the fact that d~h > 1, this proves that d s < 2d h / (1 + dh)- 
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2.3.2 Lower bound on the spectral dimension 



To prove the lower bound on the spectral dimension we exploit the following lemma derived 
in [22] in the context of generic trees and note that it can be used in a more general situation. 



Lemma 4 For all r G T one has 

Q ^ iR+ mm (29) 

The proof of this lemma can be found in [22J. It uses a decomposition of walks in fli, 
the set of walks which do not reach further than distance R from the root, and f^, the set 
of walks which do reach further. Then one can show that Q T (x) = Q^ix) + Q^ 2 (x), where 
Q^(x)<R a ndQ^(x)<^ m . 

To prove the lower bound in the theorem we observe that by existence of the Hausdorff 
dimension one has that |£> r (i?)| > R dh L(R), where L(R) varies slowly at R — > oo. Thus, by 
Lemma H] 

Q T (x) <R + x- l R- dh L~\R) (30) 
Choosing R = \x~ 1 '( 1+dh >~\ (which also optimises the inequality) proves that d s > 24/(1 + 

4)- 



2.4 Independent and identically distributed outgrowth 

If one considers random tree ensembles with a unique infinite spine, where the outgrowths 
from different vertices along the spine are independent and identically distributed (i.i.d.), 
one can make stronger statements about the dimensions of the tree ensemble. Recall, the 
definition of the hull Br as the union of the spine from r up to vertex sr and all finite trees 
attached to it. Let us furthermore denote by A 1 the union of the vertex s« and the finite 
trees attached to it. Denote by X l n the number of vertices in A 1 at a distance n from Si, e.g. 
Xq = 1. Furthermore, denote by A l n the intersection of A 1 with the ball of radius n centered 
around Sj. We have the following Theorem: 

Theorem 5 Let (r, v) be a random tree ensemble concentrated on the set of trees with a 
unique infinite spine Tf and i.i.d. outgrowths (A 1, ) i>1 on the spine. 

(i) If for a e (0, 1) and z e (0, 1), 

(z\ A \ = l-(l-z) a l(l-z), (31) 
where l(x) varies slowly at x —> + , then almost surely 

4 < - d 8 < — ^— (32) 
a 1 + a 
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provided that dh and d s exist, 
(ii) If furthermore, 



u({Xi > 0}) < 1 (33) 



where c> is a constant, then almost surely dh and d s exist and 

d h = -, d s = (34) 
a 1 + a 

2.4.1 Proof of first part of Theorem [5] 

Since \Br\ = ^2f =0 \A l \ + R is just a sum of i.i.d. random variables \A l \ (with possibly infinite 
variance), one can easily determine a bound on their sum \Br\ and with that a bound on dh 
almost surely. We begin by observing that by a Tauberian theorem [23J, ( 13TT) is equivalent 
to 

v{{\A*\ > R}) ~ R~ a L(R), (35) 

where L(R) varies slowly at R — > oo. Here, we define "~" to mean that the ratio of the two 
sides tends to one as R — > oo. By ( |35|) there exists a function L\(R) which varies slowly at 
infinity such that 

oo 

> R " L i( R )}) < °°- ( 36 ) 

i=l 

It then follows from a theorem by Feller [2U Theorem 2] that 

R 

>{{t : 1-4*1 > R ° L ^ R ) infinitely often}) = (37) 



R 

i=Q 

One then has that for //-almost all trees there exists a constant C > 1 and an Rq such that 
for R > Ro one has 

\B R \ < CR^L^R). (38) 

This proves that dh < l/a almost surely and using Theorem [TJ one has further d s < 2/(1 + a) 
almost surely. 



2.4.2 Proof of second part of Theorem [5] 

In this section we prove that there exists a constant C such that, 

u({r : \B R \ < XR^}) < Ce- x ' ai ^ lR ^ (39) 

for A small enough with l(x) from (13TT) . Assuming that (1391) holds we proceed as follows. By 
[25| Theorem 1.5], for any c£l, there exists a function L 2 (R), slowly varying as R — > oo, 
which satisfies 

um ^w-^w-'fl- 1 ^) = J (40) 

fl^oo log(i? c ) v 1 
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Therefore, by choosing A = L 2 (R) and c > 1 one finds that 

oo 

"({r :\B R \< R- L 2 {R)}) < oo. (41) 

R=l 

Using the Borel-Cantelli lemma one then has that 

v({t : \B R \ < R« L 2 (R) infinitely often}) = 0. (42) 

Thus, for //-almost all trees there exists an Rq such that for R> Rq one has 

\B R \>R^L 2 (R). (43) 

This proves that dh > 1/a almost surely. Together with the results of the previous section 
and Theorem [1] one has thus almost surely 

1 2 

d h = -, d s = — — . (44) 
a 1 + a 

The basic idea to prove ( 1391) is the following: if the finite outgrowths along the spine die 
out fast enough, i.e. if ( 133|) holds, then the ball Br and hull Br are close in size and one can 
use the latter to estimate Br. Recall that A\ is the intersection of A 1 with the ball of radius 
n centered around Sj. An essential ingredient to measure how close Br and Br are in size is 
provided by the following lemma: 

Lemma 6 For z < 1 one has for the probability generating functions 

(J A % < (z^)„ < (z\ A \ + v({X* n > 0}) (45) 

Proof The proof is inspired by ideas of [22]. Firstly, note that \A l \ > {A^l, hence the lower 
bound follows. For easy notation let us denote the event £ n = {X l n > 0}. Then 



z\ Ai \) v = [ ( z ^-z^)du+ [ (z^-z^di 
= I (z^-z^)du< I z^du 

J On J On 



< / dv = u({X % n >0}) (46) 

J £ n 



Thus the upper bound follows. 

To prove f[3"$j) we first note that since 
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it suffices to prove that 

R 



({r ■■ J2 M < XR "V < Ce- x - ai ^ lR - 1/a \ (48) 



v x k 

i=i 



Using Markov's inequality, the independence of the \A R \ and Lemma [6] we find that for any 
t E (0, 1) and A > 



<{r-J2\^ R \<XR-}) < K{r:expH^|^|)>expHA^)}) 



R 

v{{t : 1 

1=1 i=l 

R 

i=l 

< e tXRh ((e-* 1 - 411 ), + u{{X R > 0}))*. (49) 

From (I3TI) it follows that 

(e-* 1 ^ 1 ), < 1 - t a l(t) + o(t a l(t)) (50) 
where l(t) varies slowly as t — > + . Choose t = A -1 -?? -1 /" and A small enough such that 

\- a l(\- l R- 1/a ) > c (51) 

with c from ( )33|) . Apply ( |50|) and (|33|) to (1491) and get, for a suitable constant C 

i=l ^ ' 

< c'e-^"'^" 1 ^ 17 ")^. (52) 
Thus (EHD follows from (El with C = CV. 



3 Conditioned Galton- Watson trees 



3.1 The model 

In this section we apply the results of the previous section to the model of simply gener- 
ated trees. Simply generated trees can, in most cases of interest, be interpreted as a size 
conditioned Galton- Watson process. A Galton- Watson process is a discrete time branching 
process which starts from a single particle. At each time step the particles branch indepen- 
dently to k other particles with the same probability p(k), k > referred to as offspring 
probability. We denote the generating function of the offspring probabilities (or outdegrees) 
by f(z) = Y2^=o P( n ) z " ■ The value of /'(I), the mean number of offsprings, determines the 
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survival properties of the process. If /'(l) = 1 the process is said to be critical and it dies 
out with probability one. If /'(I) < 1 the process is sub-critical and dies out exponentially 
fast and if /'(l) > the process is super-critical and has a positive probability of surviving 
forever, see e.g. [26J. 

The model of simply generated trees has as parameters a sequence of non-negative 
weights (w n ) n >i referred to as branching weights and is defined by a (Gibbs) measure on 
the set of trees with n edges by 

v n {T) = Z~ 1 Yl w °(v)i TeT n (53) 
«ev(r)\{r} 

where Z n is a normalization 

Z n = ^2 Ft w "(v) ( 54 ) 

TeT n veV(T)\{r} 

often referred to as the finite volume partition function. It is useful to define the generating 
functions 

oo oo 

Z(() = J2ZnC and g(z) = J2^n + iz n (55) 

n=l n=0 

which are related by the equation 

Z(()=(g(Z(Q), (56) 

see e.g. [27J. We will denote their radii of convergence by ( and p respectively and we 
furthermore define Zq = lim^_^^ Z((). When p > the measure v n can equivalently be 
defined by a Galton- Watson branching process with offspring probabilities 

p(k) = CoWk+iZt 1 (57) 

conditioned to have n edges, see e.g. [28]. By (|56|) . the mean offspring number can be written 
as f'il) = 1 — g(Z )/Z'(( ) and thus it is clear that the process is either critical or sub- 
critical. Here we will only consider the critical case i.e. when Z'(( ) = oo in which case the 
following theorem holds. 

Theorem 7 For a critical size conditioned Galton- Watson process, the finite volume mea- 
sures v n converge weakly as n — > oo towards a measure v which is concentrated on Tc? . The 
degrees of vertices on the spine are independently distributed by 

</>(k) = ( (k - l)w k Zt 2 (58) 

and the outgrowths from the spine are independent Galton-Watson trees with offspring prob- 
abilities (5?p. 
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The critical model is usually divided into two cases depending on whether /"(l) is finite 
or infinite. For the case /"(l) < oo, Theorem [7] is originally due to Kennedy [29] and later 
to Aldous and Pitman [30J, see also [22]. To our best knowledge, the generalisation which 
includes the case f"(l) = oo was first proved in the special case w n ~ by Jonsson and 
Stefansson [2S] and later in full generality by Janson [21]. 

In the case /"(I) < oo the trees always belong to the same universality class and have 
been referred to as generic trees in the physics literature. In the infinite case there is a 
range of universality classes depending on the singular behaviour of / and this case has been 
referred to as critical non-generic. 



3.2 Dimensions 



The following theorem holds for the Hausdorff and spectral dimension in the critical case 
Theorem 8 

(i) The quenched Hausdorff dimension and spectral dimension of generic trees is almost 
surely 

4 = 2 and d s = 4/3 (59) 

respectively. 

(ii) For critical non-generic trees with w n ~ n~^L(n), where (3 G (2, 3] and L slowly varying 
at infinity, the quenched Hausdorff and spectral dimensions are almost surely 

(3-1 2(/5 — 1) 

dh = — and d s = — (60) 

h )3 -2 2/3 — 3 v ; 

respectively. 

The result on dh in (i) is proven for instance in [2D] and the result on d s was originally 
proven in [16]. The results in part (ii) were conjectured in the physics literature [321 [33] 
and later proven in the mathematical literature by Croydon and Kumagai |15j . Below we 
will show how Theorem [8] easily follows from Theorem [5j While the proof in [15] is slightly 
more general, as discussed in the introduction, the proof presented here is more intuitive 
and hopefully makes the result more accessible to physicists. 

Proof It is a standard result, see e.g. [9J, that for generic trees 

(^l) = l-c(l-C) 1/2 + o((l-C) 1/2 ) (61) 

and 

lf"(l)n + o(n), (62) 



l-f n (z) l-z 2 
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where /„ = / o • • • o /, n-times, is the n-th iterate of the generating function f(z) of the 
offspring distribution. Denoting by X]^ the size of the n-th generation of the j-th outgrowth 
from one thus has using (J62"j) that 

KK J > 0}) = 1 - /„(0) = ^-(1 + (1)), (63) 

Hence, since the outgrowths are independent, one has using Theorem [7J 

KK>o}) = i - - KW > o») 

= -(l + o(l)). (64) 
n 

The results then follow from Theorem [51 

Next, consider the case /"(l) = oo. By a Tauberian theorem [23J, u> n ~ n _ ^L(n), implies 
that one can write 

f(z) = z + (l- z )P-H 1 (l-z) (65) 

where l\ is slowly varying at zero. Let W(C) = ^(CoO/^o- Using generating function 
arguments one can deduce from Theorem [7J that 

= /'(W(z)). (66) 

Write 

w(C) = i-(i-C)^ T x(i-0- (67) 

By fl56|) and fl65|) we find that \ has to satisfy 



lim xixf^h (x^xix)) = 1 (6£ 

x-H) \ / 



which, by [251 Theorem 1.5], entails that \ is slowly varying at zero. Therefore, by (165]) . 
fl66l) . flSTj) and Tauberian theorems one can write 

(z^) = l-(l-z)Hi 2 (l- z ) (69) 

with I2 slowly varying at zero. Thus, by part (i) of Theorem [5] we have established the upper 
bounds on dh and d s . 

For the lower bounds we use part (ii) and consider the survival probability of the out- 
growths. Let X]^ be the size of the n-th generation of the j-th outgrowth from Sj. It was 
shown by Slack [311 Lemma 2] that 

KW > o}))^ 2 iiMW > 0})) = p-r^ 1 + ( 7 °) 
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where l\ is the same slowly varying function as in fl65|) . Then, since the outgrowths are 
independent, we find by fl65l) that 



which completes the proof. 

4 The attachment and grafting model 
4.1 The model 

The attachment and grafting (ag) model [10] is a recent model of randomly growing rooted 
planar trees which is a special case of a very general tree growth model, referred to as the 
vertex splitting model. The vertex splitting model was introduced in [11] as a modification 
of a combinatorial model encountered in the theory of random RNA folding [3]. The original 
motivation for studying the special case of the ag-model is that it has a so-called Markov 
branching property which makes it exactly solvable in a strong sense. It furthermore turns 
out to have a unique infinite spine which enables us to apply Theorem [5] to study its Hausdorff 
and spectral dimension. Using the first part of the theorem we establish what we believe to 
be tight upper bounds on the dimensions for the full range of parameters. The corresponding 
lower bounds require information on the extinction probability of the outgrowths and we will 
provide results on that only for a certain range of parameters. 

We give an informal description of the growth rules of the ag-model below but refer 
to [10] for a more detailed discussion. The model has two parameters a, 7 G [0,1] and a 
constraint D on the maximum degree of vertices. For easier notation we define 



Call the edges which are adjacent to vertices of degree one (besides the root) leaves and call 
the other edges internal edges. Starting from the unique tree with two edges, in each time 
step the number of edges is increased by one by randomly selecting 

(a) a vertex of degree k > 2 with relative probability i]k + 1 — 2i] — a and attaching a new 
edge to it (the k possibilities of attaching chosen uniformly at random) or 

(b) an inner edge with relative probability a and dividing it into two edges by grafting a 
vertex to it or 



KPQ > 0}) 



1 - f (1 - u({XV > 0})) 

{p - > o}))^iMPC J > o»)(i + 0(1)) 



~d + *(!))■ 



(71) 




(72) 
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(b) 

® — 




Figure 2: Growth rules of the ag-model. The root is indicated by a circled vertex. 

(c) a leaf with relative probability 1 — r\ and dividing it into two edges by grafting a vertex 
to it, 

see Fig. |2j The growth rule generates a sequence of probability measures on T which we 
denote by (V n )n>i ( n °t writing explicitly the dependency on a, 7 and D). The probability 
of a tree T G T n+1 is given by a recursion 

u n+1 (T) = "niT'MT -> T) (73) 
T'er n 

where P(X" — > T) is the probability of growing T from T" according to the above growth 
rule. 



4.2 Markov branching and convergence of the finite volume mea- 
sures 

The ag-model has a property called Markov branching, a concept originally introduced by 
Aldous [TJ. Markov branching means that for any k > 2 there is a function : — )■ M + 
such that for each finite tree To which branches at the nearest neighbour of the root to 
subtrees T 1; . . . , T k _i it holds that 

fc-i 

^Toi(To) = ftdTxl, . . . , |7W|) J] m (Ti). (74) 

i=l 

The functions qt are referred to as the first split distribution. In the ag-model, the first split 
distribution is given by [10] 



, t(ni ,...,^ l) = r(t ;" 2 ;^ r(i -" )r( "' n " r(,i -"" ) 



r(l=^T(k) rjr(n-rj) Alr(l-77)r(n i + l) 
x + (75) 



=1 4 
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where n\ H — • + rik-i = n — 1. 



In [JO] it was shown, using (TM1) and ( |75l) . that the sequence (z/ n )n>i converges weakly as 
n — > oo to a measure z/ which is concentrated on the set of infinite trees with exactly one 
infinite spine having finite outgrowths. Outgrowths from different vertices on the spine are 
independent and the probability that a vertex on the spine has degree k and that the finite 
outgrowths attached to it are T 1; . . . , T k _ 2 , is given by 

"* (r " ' ' ' ' *" 2>= r (1=*) r ( t - i)(i + m) ii r (i - ,) r ( |r,| + Ty (76) 

where m — \Ti\ + ■ ■ • + |3 r jfe_ 2 | (with the convention that ^(0) = !)• 



4.3 Dimensions 

The annealed Hausdorff dimension of the ag-model was calculated in [10] using a generating 
function argument. It was shown that du = 1/ot when D is finite and in the case D = oo 
and a > 1 — 7. Note that it is meaningless to consider the annealed Hausdorff dimension 
when D = 00 and a < 1 — 7 since then it follows from (1761) that the expected degree of a 
vertex on the spine is infinite and thus dn = 00. 

Despite this, one might still conjecture that a.s. d^ = I /a. This was confirmed in 
[TO] when 7 = in which case the outgrowths from the spine are single edges and it was 
furthermore shown that in this case d s = 2dh/ (1 + dh) = 2/(1 + a). Below, we extend these 
results to a wider range of parameters. 

Theorem 9 1. For the ag-ensemble (T, v) it holds for any a, 7 and D, that almost surely 

1 2 

dh < — and d s < . (77) 

a 1 + a 

2. Furthermore, for D = 3, h < a < 1 and for D = 00, 7 = 0, we have almost surely 



1 2 

dh = — and d s = . (78) 

a 1 + a 

Proof As before, we denote the outgrowths from vertex Sj on the spine by A 1 . The upper 
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bounds in f J77|) follow from the simple observation that 



oo oo 



& A X = E E E E 

fc=2 m=k-2 mH hn fe _ 2 =m Tier ni ,...,T fc _ 2 er„ fe _ 2 







1-a 



°° ar fc - 2 + 

V / r ^ [1 - (1 - X) 

^ r (if) r(* - 1) 



fc=2 

1 - (1 -2)° 

Z 



(79) 



The second equality in (1791) is obtained by noting that the innermost sums of u ni (Ti) over 
Tj, i — 1, . . . ,k — 2, yield 1 and then using standard gamma function identities [35J 



EaTin — a) , , „ 

rd-W+i) ' =1 - (1 - s) - (80) 

n=l v ' v ' 

Finally, an application of the first part of Theorem |5] yields (I77|) . 

To obtain the equalities in (175)1 we study the height distribution of the outgrowths from 
the spine and apply the second part of Theorem [5j We start by noting that the case when 
D = oo and 7 = 0, which was already proved in [10J, follows immediately from Theorem |5] 
since in this case the outgrowths from the spine all have height 1. We now focus on the case 
D = 3. As before, denote by X % n the size of the n-th generation of A 1 . The result follows 
from 

Lemma 10 When D = 3 and | < a < 1 it holds that 

KK>o})< + (si) 

n 

We devote the following subsection to the proof of this lemma. 



4.4 Proof of Lemma (TH 



We start the proof with a general approach and reduce to the specific parameters stated in 
Lemma [10] only in the end when needed. In this way we keep the problem of extending the 
results to all parameters clearly approachable. 

For a finite tree T, let h(T) be its height, i.e. the maximum distance from the root to 
any vertex. We make the following definitions for more compact notation 

c k = T ( k ~\ +1 ^ and CnB =-Vnn-vW({h(T)^R}) 

r(!=s)r(ifc) r(i-tj)r(n + i) 1 ; 
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with the convention that Vo({h(T) < R}) = 1. The main tool that we will use in the proof 
is the following generating function 



H R (C) = J2 Cu,rC- 



n=0 



Proposition 11 The probability that A % is extinct at level R is given by 



u({X R = 0}) = a / (H R (()) — d(. 
Jo 

Proof Using the distribution of the outgrowths given in (ITB"]) one can write 
u({X R = 0}) = u({h(A*) < R}) 



a 



34) 



oo oo k—2 

E( fc - E E U(- G ^ ni d ^ ( 85 ) 

. fc=2 m=k— 2rtiH hrifc-2= m > n i>l i=l 

The sum over m is first performed giving (1 — H R (Q) 2 and the sum over k is then performed 
using standard gamma function identities similar as above giving <\84L 



Proposition [TT1 shows that H R (() contains all the information needed to calculate the extinc- 
tion probability of A 1 . Next, using the Markov branching property, we will derive a recursion 
equation for H R (Q which enables us to extract enough information. 

Proposition 12 For R > 1, it holds that 



OL-\-7} — 1 



H R (C) = l + aH R -x(QJ {H R ^{y))—dy-{a + ri) / (H R ^(y))vdy (86) 
with the convention that H (() = 1. 



Proof Using the Markov branching property and (1751) we write 

oo k— 1 

-nC n . R = ^c fc Y Yi{-C ni ,R-i) 

k=2 niH \-n k _ 1 =n— l,n»>l i=l 

/ fc-1 

x I 1 — ?7 — a + a 



n — n~ 



57) 
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Multiply by £ n and sum from n — 2, . . . oo to get 

oo oo k— 1 

-too +^ = E c *E E n(-^*-iK ni 

fc=2 ?i=fc «iH hifc_i=n— l,r»i>l i=l 

/ fc-1 

x I 1 — T) — a + a 



3=1 3 
x / \ fc— 1 



' ■ \ fc 2 



k=2 \ n=l , 

(, OO \ OO / iaj \ 

- ]T a^-iC I jf - x ) c * ( " E ^-i*/ 71 j dy. (88) 

Using the definition of H R and performing the sum over k yields 

H' R (C) = -v (Hn^C))^ + «#W0 / (H R ^(y))^ dy. (89) 

Jo 

Integrating this equation and using integration by parts on the second term on the right 
hand side finally yields the result (|86|) . 

We now reduce to the case D = 3 in which case rj = a — 1 < 0. The recursion (1861) is then 
greatly simplified. We define F R (() = a f£ H R (y)dy and then y{{X R = 0}) = F R (1). Next 
we integrate ( 1861) to get 

^(C) = «C + -i^-i(C) 2 - — / / {F^ix)) 2 dxdy. (90) 

2 a Jo Jo 

One can see directly from the definition of H R fl83|) that when 77 < 0, Fr and all its derivatives 
are increasing in R. One can furthermore deduce from (I8"3"|) that 

Foo(C):= limF*(C) = 1- (1-Cr- (91) 

i?->-oo 

Therefore, under the assumption that | < a < 1 we may insert -Foo(C) into the integral in 
(|90|) to get the inequality 

^(0>^ + ^(^-i(C)) 2 -^(l-C) 2a . (92) 
Evaluated at ( = 1 this can be written as 

i^(l) > /(^-i(l)) (93) 

with f(x) = \ + \x 2 . We can interpret / as a generating function of the offspring probabilities 
p(0) = p(2) = 1/2 of a critical Galton- Watson process, cf. Section [31 Let f R be the i?-th 
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iterate of / i.e. fa — f o ■ ■ ■ o /, R times. Now, / is an increasing function and therefore, by 
repeatedly applying fl93|) one gets 

p({X R = 0}) = F R (1) > f R - 2 (F (l)) > / B _ 2 (l/2) (94) 

where in the last step we used that Fq(1) = a > 1/2. It then follows from ( J62l) that 
/r(1/2) = 1 — -|(1 + o(l)). This concludes the proof of Lemma [TU1 

5 Discussion 

We developed relatively simple methods, relying on generating function arguments, for cal- 
culating the Hausdorff and spectral dimension of trees with a unique infinite spine resulting 
in Theorems [1] and [51 where the latter applies when the outgrowths along the spine are inde- 
pendent and identically distributed. These methods were applied to two models of random 
trees of very different nature, demonstrating the versatility of our results. 

The first application of Theorem [5] concerns non-generic critical trees which are a special 
case of simply generated trees. The values of the Hausdorff and spectral dimension of the 
non-generic, critical, size conditioned Galton-Watson trees were conjectured by mathemat- 
ical physicists 15 years ago (321 [33] but only recently proved by mathematicians [15]. We 
used the opportunity here to communicate these results to the physics community as well 
as providing a simple alternative proof. 

We would like to point out that these results complete the study of the dimensionality 
of the different phases of simply generated trees. Theorems [7] and [8] concern critical Galton- 
Watson processes. As was mentioned in Section [3X the simply generated trees can also 
correspond to sub-critical Galton-Watson processes. In the sub-critical case, the finite 
volume measures converge weakly to a measure concentrated on the set of trees with a 
finite spine ending with a vertex of infinite degree. The length of the spine has a geometric 
distribution and the outgrowths from the spine are finite, independent, sub-critical Galton- 
Watson processes. This convergence theorem was originally proved in [28] in the case w n ~ 
n" 13 and later in full generality in [31 J . 

Due to the presence of a vertex of infinite degree it clearly follows that almost surely 
d s = oo since a random walker will eventually hit the infinite degree vertex and similarly 
dh — d>H — oo. However, it was shown in [2S] that when w n ~ n _/3 , the annealed spectral 
dimension is finite and ds = 2((3 — 1). This is due to the fact that the fluctuations of the 
outgrowths from the spine serve to slow the random walker down on the way to the infinite 
degree vertex. Note that in this case it does not hold that ds = 2<i#/(l + du) as in the other 
phases due to the absence of the infinite spine. 

We summarize the above discussion in the phase diagram in Figure [3] where we consider 
the case when w\ and j3 are the free parameters and w n ~ n _/3 . This choice of parameters 
allows us to access all phases and explore the full range of dimensions in each phase, see [2B] 
for a more detailed explanation. 
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Figure 3: A phase diagram for the size conditioned Galton- Watson trees, having free pa- 
rameters Wi and 0, where w n ~ n~P . A critical line separates the generic phase from the 
non-generic phase. For the values 2 < /3 < 3 on the critical line (grey line) the offspring prob- 
ability distribution has infinite variance and therefore this case corresponds to non-generic, 
critical trees. 

The second application of Theorem [5] concerns the attachment and grafting model which 
is a special case of the vertex splitting model. We give a novel proof that for the parameter 
range D = 3, § < a < 1 and for D = oo, 7 = 0, one has almost surely dh = 1/ot and 
d s — 2/(1 + a). This proves part of a previous conjecture [TO] . 

It remains an open problem to prove that the quenched Hausdorff and spectral dimension 
of the ag-model are almost surely d^ = 1/a and d s = 2/(1 + a) for the full range of the 
parameters a, 7 and D. The only missing ingredient in the proof is to generalize Lemma [TU] 
to the full range of parameters. We formulate this as the following: 

Conjecture For the ag-ensemble (r, v) it holds for any a, 7 and D, that there exists a 
constant c > such that 

KK>0})< -(l + o(l)). (95) 
n 

Furthermore, we expect that one even has the stronger result ^({X^ > 0}) ~ I /(cm). This 
conjecture could possibly be proven by analysing the recursion relation fl86|) to find the 
critical behaviour of the generating function Hr(£) when £ — > 1 and R — > 00. We hope to 
return to this proof in the near future. 

At this point we would also like to note that all of the results derived for the ag-model 
can straightforwardly be extended to Ford's a-model [2] and its generalisation, the 07- 
model [36]. We leave a more detailed discussion for future work after the proof of the above 
conjecture. 
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There are also different applications of the presented work in the field of quantum grav- 
ity. Firstly, the results derived in the article for critical non-generic trees are useful to 
understand in more detail the branched polymer phase of Euclidean quantum gravity and, 
in particular, in the presence of matter. Furthermore, the critical non-generic trees are 
described by branching processes whose population size has polynomial growth which is po- 
tentially relevant to the analysis of recently introduced multigraph ensembles as models for 
four-dimensional causal or Lorentizan quantum gravity [37J [38] . 
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